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Abstract 

We consider models of accelerated cosmological expansion described by the Einstein 
equations coupled to a nonlinear scalar field with a suitable exponential potential. We 
show that homogeneous and isotropic solutions are stable under small nonlinear perturba- 
tions without any symmetry assumptions. Our proof is based on results on the nonlinear 
stability of de Sitter spacetime and Kaluza-Klein reduction techniques. 
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1 Introduction 

At present the subject of accelerated cosmological expansion is of great astrophysical interest. 
Many candidates have been suggested for the cause of the acceleration, known under the name 
of dark energy. The greater part of the literature on this concerns homogeneous (and even 
isotropic) models or their linearized perturbations. This is often enough to make contact 
with observations. Nevertheless, since many of the phenomena of interest are linked to 
inhomogeneities, it is desirable to develop an understanding of the full nonlinear dynamics 
for initial data which are as general as possible. One approach, which is pursued in this 
paper, is to try to prove general mathematical theorems. A recent review of this approach 



A well-known feature of models for accelerated cosmological expansion is that they exhibit 
attractor solutions which are homogeneous and isotropic. The simplest model is a positive 
cosmological constant in which case the attractor is the de Sitter solution. In that case a 
theorem on stability of this solution under small nonlinear perturbations has been proved [2] • 
It concerns the vacuum equations with positive cosmological constant. No symmetry as- 
sumption is required. Generalizations where matter such as a perfect fluid or kinetic theory 
is added are not available. Under the restriction of plane symmetry an analogous result with 
collisionless matter was proved in Jl] . There are also no generalizations of the result of 5 
in the literature to other models of accelerated expansion such as nonlinear scalar fields. 
This paper proves a generalization of this type for a nonlinear scalar field with exponential 
potential. Because it was originally applied in models of the early universe this is associated 
with the name power-law inflation. The method of proof restricts the exponents allowed to 
a discrete set. 

The result of [5| makes essential use of the conformal properties of the Einstein equations in 
four dimensions. Anderson has extended this analysis to to any even spacetime dimension. 
Here we combine the results of with a Kaluza-Klein reduction which relates power-law 
inflation in four dimensions with a suitable exponent in the potential to vacuum spacetimes 
with cosmological constant in higher dimensions. It is proved that certain homogeneous and 
isotropic solutions are nonlinearly stable. 

Consider a spacetime (M,g,<p) that satisfies the Einstein equations 



where V{(p) is the potential of the scalar field. Using that the spacetime is four-dimensional, 
Eq. (Q) and (j2j) can be condensed into 



is [H]. 





with nonlinear scalar field matter 




(2) 





The Bianchi identity implies the equation of motion for the scalar field 



Dip = V'(<p) . 



(3b) 
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Power-law inflation refers to the case 

V(ip) = V exp [— kA if] (4) 

with constants Vq and A, where A < y2. Models with a potential of this type are the subject 
of the following. 

The most elementary power-law inflation models are the homogeneous and isotropic mod- 
els [J] , of which the simplest is 

2+ — 

g = -dp + (jpj d P + i Stj di l dS: j , (5) 

where d is related to the exponent A in the potential via A 2 = 2d(2 + e£) . The main theorem 
of this paper establishes nonlinear stability of this model. 

Theorem 1. Consider the exponential potential 



V(<p) = Vq exp f-M/2y - — - ^pj with d G N, d even . (6) 

On T 3 consider smooth Cauchy initial data for the Einstein scalar field equations © with 
potential ©. Let the initial data be close in the C°° -topology to the homogeneous and isotropic 
initial data of the model (J5J). Then the initial data evolves into a power-law inflation model 
(M, g, ip) which is geodesically complete to the future, and there exists a coordinate system 
(t,x l ) that is global to the future, such that g takes the form g = —a 2 dt 2 + gijdx l dx^ , where 
gij and ip admit the asymptotic expansions 



l d + 2 

m>0 m>0 



w = p+4/d E 4 m) l ~ 2mld > ( p = «- 1 ^\l^r log * + S ^ i ~ 2m/d > ( ? ) 



and a = Ylm>0 a ( m ) ^ 2m / d . Thus, the homogeneous and isotropic power-law inflation model 
(jSJ) is nonlinearly stable. 



A slightly different formulation of the theorem uses the concept of asymptotic simplicity. We 
call the spacetime (M,g,(f) asymptotically simple (in the future), when g is conformal to 
a metric g = Q 2+d g, and if = if + c?k _1 A _1 log S7, where g, ip, and the positive function fi 
can be smoothly extended through the hypersurface f2 = 0, which we denote as the 
conformal boundary of M. The formula for the conformal transformation of the curvature 
scalar implies g^V^QV u n = -2kT 2 Vq(2 + d) _1 (3 + d)~ l exp(-KA0) on J^+, hence is 
spacelike. Note in this context that our definition for asymptotic simplicity differs from the 
standard definition which applies for the vacuum case and for the case when matter can be 
neglected in an appropriate sense in the neighborhood of the conformal boundary. 

Theorem n states that initial data close to homogeneous and isotropic data evolves into an 
asymptotically simple solution. We also have 

Theorem 2. The initial value problem for the Einstein equations with scalar field matter 
and potential (jfi]). where initial data is prescribed at conformal infinity J"^, is well-posed. 
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The resulting spacetime (M,g,cp) is asymptotically simple with conformal boundary ^ + and 
the metric and the scalar field admit the asymptotic expansions Ffy, where the coefficients are 
uniquely determined by the initial data. 

The proof of the theorems is based on the fact that power-law inflation models (M, g, <p) with 
potential © are in one-to-one correspondence with a certain type of d-dimensional vacuum 
solutions {M,g) of the Einstein equations with positive cosmological constant. This is proved 
in Section |2 In Section |3] we give a brief overview of the asymptotic behavior of solutions 
(M, g); in particular we recapitulate existence and nonlinear stability of asymptotically simple 
solutions. On the basis of these results, the theorems are proved in Section 

Theorem ^ formulates the asymptotic behavior in a certain coordinate system that is not 
Gaussian in general. In Appendix El we briefly discuss asymptotic expansions in Gaussian 
coordinates; in particular we show that this choice of coordinates introduces logarithmic 
terms into the expansions. In Appendix El we investigate asymptotic expansions for general 
exponential potentials on the level of formal power series. We demonstrate that logarithmic 
terms that appear in the series can often be removed by a suitable choice of (non-Gaussian) 
coordinates; this complements previous studies of formal expansions of this type [Sj. 

2 Reduction 

In this section we employ the Kaluza-Klein reduction method to find a one-to-one relationship 
between solutions of the Einstein equations with positive cosmological constant A and certain 
power-law inflation models. The formulation of Kaluza-Klein theory used and the notation 
is primarily based on j^j. 

Consider a principal fiber bundle G — > M —> M , where the base space M is a 4-dimensional 
differentiable manifold, and G a <i-dimensional Lie group (which we will eventually assume 
to be abelian). On M, let g be a Lorentzian metric such that g is invariant under the right 
action of G and vertical vectors are not null w.r.t. g. The metric g induces a Lorentzian 
metric g on M, a metric £ on each fiber, where £ is invariant under the action of G, and a 
connection on M (in the form of a horizontal bundle). In the so-called polarized case the 
horizontal distribution is assumed to be involutive and in an adapted local trivialization over 
a chart neighborhood of M with coordinates {x^ \ fj, = . . . 3} the metric g can be written as 

9 = 9AB e A e B = 9fiu dx^dx" + £ mn 9 m 9 n , (8) 

where the 9 m constitute a basis of right-invariant 1-forms on G, and gnu, £,mn depend only on 
the coordinates {x 5 }. Greek indices run from to 3, latin indices m, n, etc. assume values 
1,2, ... ,d. Capital letters A, B, etc. run over the combined range: the components Tab of a 
tensor thus comprise T^ v , T^ n , T mu , T mn . 

We set 

:= ^/dete so that Un = e 2(j,/d ( mn , (9) 

where det £ = 1. Like £ TOn , in the given trivialization, the field <j) can be regarded as a scalar 
field on M. 
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To compute the Ricci tensor Rab of the metric g we recall the general Kaluza-Klein formulas 
from [3] and use V M <^> = ^ mn V M ^ mn , where £J nn is the inverse of £ mn . We obtain 



Rmn — Rmn 2~£,^^ a^mp^ ^ n Q 2 ^ Q ^ ^^q) £r. 



(10a) 
(10b) 



where R^ u is the Ricci tensor of g^ iu and R mn the Ricci curvature of the fibers. Both R ri 
and the components R^ n of the Ricci tensor Rab vanish if the Lie group G is abelian. 

On M we introduce the conformally rescaled metric 



9fiu '■ — ^ 9[iu 



(11) 



and denote the Ricci curvature of g^ u by R a y Employing V a V a + V Q 0V Q = V a V a , 
where V a = g a/3 V/3, Eq. (|T0|) becomes 



1 



1 1 



i?^ = fl^ + - ( V 5 V> ) ^ - ( - + - ) V^V U (P + -V^nmVuC 



1. 



/? — /? — -p 9 



2 ,~ 



(V a V<V)£ mri + e 2 ^ d [V a V a Un - ( pq V a ( mp V a C 



nq 



(12a) 
(12b) 



Contracting (|12bl) with £ mn leads to 

r^mn = H mn Rmn ~ ^ V a V> , (13) 

where again £ mn i? mn = in the case of an abelian Lie group G. Eq. 1)12(1 simplifies when £ mn 
is independent of x s , i.e., when V M Cmn = 0. 

Assume that (M,g) satisfies the Einstein vacuum equations with cosmological constant A, 
i.e., 



Rab 



d + 2 



A 9AB ■ 



(14) 



Suppose further that the Lie group G is abelian and that V^Cmn = 0. From (|T3|) we obtain 

2d 



d + 2 



Ae~ 



Eq. (|12a|) then leads to 



Define 



n t ,r - Ac-'.,),,,, + ( ~ + ~ ) \>v, 



-1 



(/3 = K 

then (|T5)) and (|TB|) become 

V«V> = y'(9?) and R^ = K 2 V(ip) g^ u + k 2 V^V u <p 

where 

V{ip) := Ak" 2 exp 



-k\/2 



d + 2 



(15) 

(16) 
(17) 

(18) 
(19) 
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By comparing (]18[) with Q3aj) and ()3b|) we conclude that (M, 5, is a solution of the Einstein 
equations with nonlinear scalar field, where the potential is an exponential function, i.e., 
(M , g, ip) is a power-law inflation model. The exponent in the potential (|19|) is 



\ = \ d := ^2y^ < V2, (20) 

cf. ©. 

Conversely, given a solution (M,g,ip) representing power-law inflation with exponent A = 
\<i = (2d) 1 ^ 2 (d + 2) -1 / 2 for some d G N, we are able to construct a (4 + ci) -dimensional 
solution (M , <?) of the Einstein vacuum equations with positive cosmological constant. We 
take for G an abelian Lie group, which ensures R^ n = 0, and we set d^C mn = 0; e.g., we use 
M = M x G where G = M. d is endowed with ( mn = 5 mn . The relations (|18|) together with 
the form of the potential imply Rab = 2/(d + 2) A cjab with A = k 2 Vq. 

We conclude by giving a schematic overview of the one-to-one correspondence of solutions, 
which has been established: 

(M = M xR d ,g) (M,g,tp) 

V - = n 2 T^[V(tp)] , Dip = V'(ip) 

Rab - \Rqab + ^9ab = 



(21) 



V((p) = Ak 2 exp -K\/2J-^ip 



g = g^ dx^dx" + e 2 ^ d S mn dy m dy n g^ = e* , <p = k' 1 ■ 



2d 



For the de Sitter solution in (4 + d) dimensions we can write 



g = -( dx y + e 2tlx (dx L y + (dx 2 y + (dx 6 y + (dy l y + . . . + {d y d y , (22) 

where H~ 2 = (d + 2)(d + 3)/(2A). From J2D 

we infer that (f) = dHx°, and g becomes 
g^dx^dx" = e dHx ° (-(dx ) 2 + e 2Hx ° [(dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ]) . (23) 
By introducing new coordinates x^ through dx° = exp(<iffx /2) dx° and x l = x % we obtain 

g^dxV = -(dx ) 2 + +d (x f + ~ d Sij d^dx j , (24) 

i.e., a flat Robertson- Walker model for power-law inflation as in ijH]). 



3 Asymptotic series 

Consider the Einstein vacuum equations with cosmological constant A in n + 1 dimensions, 
n > 3, n odd. The n + 1 decomposition of the equations consists of the constraint equations 
and the evolution equations 

2A 

dt g ab = -2g ac k c b , d t k% = R\ + (txk)k% -5% , (25) 

n — 1 
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where we have used a vanishing shift vector and a lapse function set equal to one. In jllj it 
was proved that the equations admit power series of the following type as formal solutions, 

9a b = <? m (€ + e~ 2Ht € + e~ 3m + • • •) (26a) 
k% = -HS% + £ fcV) e ~ mHt > ( 26b ) 

m>2 



where H = y / 2A/[n(n — 1)]. The coefficients k\ m) = o" a 6(m) + n 1 trfe (m) 5 a b are obtained 
recursively through the relations 

m— 2 m—2 

[n - m]Ha a Hm) = <r a &(p) trA: (m _ p) + tf R a b (m) , [2n - m]Itrfe (m) = ^ trfc (p) trfc (in _ p) + i? (m) , 

p=2 p=2 

(27a) 

for m > 2, m ^ n, m ^ 2n, which follows from (|25|) . and through 

m-2 

2(n - l)mrA (a) = R (ra) + ^ [-* a 6(p A(-P) + trfc (p) trfc (m _ p) ] , (27b) 

p=2 

for m = 2n, which follows from the Hamiltonian constraint, k , , , vanishes for all odd m < n. 



The evolution equation dt g a b = — 2g ac /c c b implies that the coefficients g^T' are determined by 

a 

6(1)' 



the coefficients fc a 6(1) , I = 0. . . m, recursively; in particular = for all odd m < n. The 



remaining unspecified coefficients and encode the free data, 

€=A ab , 9% = B ab , (28) 

where A ab is a Riemannian metric, -B a (, a symmetric tensor that satisfies A ab B a b = and 
^ a B a b = 0, where A ab is the inverse of A a b and V a refers to A a b. 

Consider now a spacetime (M,g) that is asymptotically simple and de Sitter (in the future), 
see, e.g., jSj. By definition, g is then conformal to a metric g = Sl 2 g, where g and the positive 
function 0, can be smoothly extended through the hypersurface SI = 0, which is often denoted 
as the conformal boundary of M. Since g^ v V ^VLV ' v Vl = —2An~ 1 (n — 1) _1 = — H 2 on 
which follows from the conformal transformation of the curvature scalar, the metric g 
takes the form 

g = -H- 2 a 2 dQ 2 + g ab dz a dz b , (29) 

when f2 is used as the time coordinate and the z a are spatial coordinates that are constant 
along the curves orthogonal to slices f2 = const. The function a depends on Q and z a ; 
a = 1 on J? + . Letting t = exp(-Hft) the physical metric becomes g = —a 2 di 2 + g a \,dz a dz b 
with g ab = exp(2Ht)g a f,- In [111 Sec. 4], in dimension n = 3, it is shown that this relation 
together with an analogous set of fall-off conditions for a and k a b implies that Gauss coor- 
dinates can be introduced in which the metric g a h and the extrinsic curvature k h exhibit an 
asymptotic expansion of the form (|26[) . It is straightforward to apply the proof in to all 
odd dimensions. 

The initial value problem for the Einstein equations with positive cosmological constant, 
where initial data is prescribed at conformal infinity is well-posed; this has been shown 
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in [3] in the case n = 3. In particular, given an arbitrary Riemannian metric A ab on a 
(compact) manifold and a symmetric tensor B ab that is tracefree and divergence- free, 
then there exists a unique future asymptotically simple solution of Einstein's equations whose 
conformal boundary is J^ + . Global non-linear stability of asymptotic simplicity has been 
proved in [§]. Hence, the evolution of initial data sufficiently close to de Sitter data yields 
a spacetime that is globally close to de Sitter and in particular asymptotically simple in the 
past and in the future. For our purposes it is most relevant that these statements have been 
generalized recently to arbitrary odd (spatial) dimensions n, in particular, even-dimensional 
de Sitter spacetime is (globally) non-linearly stable, see pQ. 

We conclude that any initial data close to de Sitter evolves into an asymptotically simple 
solution having an asymptotic expansion of the form given in (|2(i[) together with ([28JI , where 
A a b and B ab correspond to the conformal initial data set. 



4 Reduction of asymptotic series 

Consider the (4 + d)-dimensional manifold M = M x M. d ; let d be even. On M consider 
solutions of the Einstein vacuum equations with cosmological constant A of the type 1)21(1. 



g = -dt 2 + g ij dx i dx j + e 2(t>/d 5 mn dy m dy n . (30) 



g ab dz a dz 



a r ]~b 



If (M,g) is asymptotically simple, as is guaranteed for solutions sufficiently close to de Sitter, 
g a b exhibits the asymptotics (|26|) , It follows that 

9V = e 2Ht \9<$ + £ gff e- mHt ) , <t> = dHt + <t> (0) + £ * (m) e~ mHt . (31) 



m>2 / m>2 



From the 3 + d split of k a b , 

k% ^®dz b = k\ ±. ® cfc* + x5Z JL ® d y- , (32) 
we obtain in an analogous manner 

k i j = -HS i j + ^ f k i j(m) e- mHt , x=-H+Y J X^z- mHt ■ (33) 

m>2 m>2 

The evolution equation dt g a b = —^9ack c b , cf. (|25[). reduces to 

<9 4 = - 2 9ijk\ , <9 t = -dx , (34) 



hence <^™ 5 (m > 2) is determined recursively from gf) and i = 2...m, and x, 



(m) 



8 



J.M. Heinzle and A.D. Rendall Power-law Inflation in Spacetimes without Symmetry 



mHd V (m) . Reduction of the recursive algebraic system (|2*7|) yields 



m-2 



[3 + d - m]flV i(m) = "P) (m) + I tr V P ) + dx (m . p) ] (35a) 

p=2 
m-2 

[(6 + d) - m]Ftr/s (m) + 3Hdx (m) = P {m) + [trfc (p) + dx (p) ] trfc (m _ p) (35b) 

p=2 
m-2 

dHtik (m) + [(3 + 2d) - m]Hdx (m) = p (m) + ^ [trfc (p) + dx (p) ] dx (m _ p) , (35c) 

p=2 

P*, = ^ - ^VVV,> - V l V^ = 2 P j (m) e- mHt (36) 



m>2 



and p = — V l Vj<^> — V'^Vj^ with an analogous expansion. (For the reduction it is useful 
to employ the warped product structure of the metric, see, e.g., [4-.) The m th coefficient 
P 1 a (m) is determined by the coefficients gfj and cf) m , with / = ... (m — 2) ; p {ui) by (j) m , 
I = 0. . . (in — 2). The determinant of the coefficient matrix of the l.h.s. of ()35bl35cj) is 
(m — [d + 3])(m — 2[d + 3]). The system (|33|) thus determines ^ l j {in) , trfc (m) , x (m) recursively 
except for m = d + 3, m = 2(d + 3). In the case m = 2(d + 3) the system ()35bl35c[) is 
complemented by the equation 

m—2 2 m—2 

2H(2 + d) (trfc (m) + dx (m) ) = - ^ ^j ( p)<^ (m . p) + ^ trfc (p) trfc (m . p) + 

p=2 p=2 
, ^ m—2 m—2 

H -j— ^ dx( p) dx {m _ p) + 2 ^ dx (p) trk (m _ p) + P (m) + p (m) , (35d) 

p=2 p=2 

which is the reduced version of the constraint equation ()27b|) . 

The th and the (3+d) th coefficients are undetermined by the algebraic system, they represent 
the free data. When we decompose A ab , B ab , cf. (J2EJ), according to 

A ab dz a dz b = A ij dx i dx j + e 2A 5 mn dy m dy n , B ab dz a dz b = B ij dx i dx j + B5 mn dy m dy n (37) 

we obtain 

g^=A lJ: <t> (0) =dA, g% +d) =B tJ , (e 2 ^) (3+d) = B . (38) 
Hereby the data must satisfy the following conditions: 

A'-'llij + dBe~ 2A = , V'Bij - dBe' 2A VjA + dBijA ik V k A = , (39) 
where A %3 is the inverse of Aij and Vj refers to A^. 
We now make use of the relation (121[) to prove Theorems Q & [2 

An asymptotically simple solution (M, g) of the type (|3*0|) uniquely corresponds to an asymp- 
totically simple solution (M,g,(p) representing power-law inflation, 



~ g = ^> (-dt 2 + g lj dx i dx j ) , <p = K 1 J~^f<f> 



(40) 
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Thus, the asymptotic behavior of and ip is uniquely determined by studying the (reduction 
of the) asymptotic expansions of g. From the above analysis we obtain that the asymptotic 
behavior of g^ u and ip is given by the asymptotic series 1)3 ljl of g^ and <j). The coefficients in 
these series are determined via 1)34)1 through the coefficients and which are in turn 
determined recursively by ()35)) , The remaining free data is specified as in 1)38)) . 



Introducing a new time coordinate t through dt 
metric ()4U)I becomes 



exp(dHt/2) dt and setting x % 



where 



-dm) d p + e {4>~dHt) 



-or 



dH 



72- 



"3 /l. 



'J 



dx % dx 3 , 



the 



(41) 



a = exp 



m>2 



(m) 



m>2 



IT 



-2m /d 



-2m/d 



(42a) 
(42b) 



To show Theorem O consider on the three-dimensional manifold a Riemannian metric 
Aij, a symmetric tensor Bij, and fields A, B, that satisfy the condition ()39)) . Defining j4 a b, 
-B a fc according to ()37)l results in Cauchy data at conformal infinity = x IR d for the 
Einstein vacuum equations with positive cosmological constant in (4 + d) dimensions. The 
well-posedness of the corresponding Cauchy problem has been established in pQ. Reduction 
of this result yields Theorem |^1 

It is straightforward to show from 1)35)1 that 1)41)1 together with 1)42)1 coincides with the ho- 
mogeneous and isotropic solution Q when A = 1, Aij = and (e 2< ^ //a! )(3 + d) = B = 0, 
<^i +d) = 5jj = 0. Since the solution © uniquely corresponds to the (4 + (i)-dimensional 
de Sitter solution, nonlinear stability of the latter reduces to nonlinear stability of the for- 
mer, which shows Theorem ^ Equivalently, Theorem |21 can be applied directly to obtain 
the result. Hence, the evolution of initial data sufficiently close to data characterizing the 
power-law inflation Robertson- Walker model yields a spacetime that is globally close to that 
model and the spacetime possesses a metric g and a scalar field <p of the form (|4Uj). which 
exhibit the asymptotic expansion ()3T)) in a future end of the spacetime. 



5 Conclusions 



In this paper it has been shown that in certain models of accelerated cosmological expansion 
homogeneous and isotropic solutions are stable under small nonlinear perturbations with- 
out any symmetry assumptions. These results concern the Einstein equations coupled to a 
nonlinear scalar field with a suitable exponential potential. They show that some known 
results for spacetimes with positive cosmological constant generalize to a situation where the 
acceleration of a cosmological model is due to the effect of a nonlinear scalar field. 
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For cosmological applications it would be desirable to incorporate a description of ordinary 
matter (galaxies and dark matter) into the models. It is expected that the source of the 
cosmological acceleration (dark energy, the cosmological constant or the scalar field) will 
dominate the dynamics at late times, but this should be proved rather than assumed. In 
this paper we were not able to include ordinary matter but note that this has not yet even 
been done for a perfect fluid or collisionless matter in the case of a cosmological constant. It 
seems that in order to do this methods will be needed which are more direct than those using 
conformal invariance properties. 

Another direction in which the results should be extended is to nonlinear scalar fields with 
more general potentials. The case of an exponential potential with a general value of the 
exponent is discussed on the level of formal power series in Appendix B of this paper. The 
observation that a judicious choice of time coordinate can simplify the asymptotic expansions 
may be useful for later work using other methods. A similar discussion for a potential 
with a strictly positive lower bound is given in For wider classes of potentials the only 
mathematical theorems concern spatially homogeneous spacetimes of Bianchi types I- VIII, 
including normal matter jSJ I1UI 1 1 3 j . 

Questions related to cosmic acceleration and dark energy play a key role in modern cosmology. 
They deserve the attention of researchers in mathematical relativity and we hope that this 
paper will contribute to the development of this area of mathematical physics. 



A Asymptotics in Gaussian coordinates 

The coordinates in which the asymptotic expansions have been given above are not Gaussian; 
in this section we investigate the asymptotic expansions in Gaussian coordinates. We begin 
by showing that the spacetime admits Gauss coordinates that are global to the future. 

The metric and extrinsic curvature functions satisfy the following estimates: 

\9ij\ < CP +4/d , \g ij \ < Cr 2 - A ' d , |fj. fc | < C (A.la) 
| Q _ e «U/2| < ci -A/d ; |^ a | < crf-l-4/d ? < q ( A _ lb ) 

lefy < Ci~ l -^ d , |trjfe + (3 + 6/d)e~ dA/2 t' 1 ] < Ct~ x ~^ d (A.lc) 

Lemma 3. Consider a metric g iiv dx il dx v = —a 2 dt 2 + gijdx l dx^, cf. 141}) , which is given on a 
time interval [T, oo), and assume that there exists CsR such that the estimates \A . 1\) hold. 
Consider a hypersurface t = to and Gaussian coordinates based on that hyper surf ace. If to is 
sufficiently large, then the Gaussian coordinates extend globally to the future. 

Proof. Consider an affinely parametrized geodesic 7(7") that is orthogonal to the hypersurface 
t = to. By a slight abuse of notation we write (t,x)(r) for 7(r); 

t(T ) = t , x(t ) = x'q , —(r ) = a~ x (t ) , — (r ) = . (A.2) 
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It is a solution of the geodesic equations 



d 2 t _-■ / dt\ 2 dx 1 dt i ~ dx 1 dx 7 

— 2+ a- d~ t a f - ) + 2a- d t a—- + a % — — 



d 2 x l „■ f dt\ 2 2 r dx l (if , dx J df -,• dx 1 dx 7 ,. , . 

+ a8 a UJ " 5 Qtrf ^S " W J*5F + p * * * " ' (A ' 3b) 

Let e > be sufficiently small in comparison to a _1 (ro) and E > 0, let ( G (1 + 2/d, 1 + 4/d), 
and consider the maximal interval [to,t) such that 

\^-a~ l \<e, $&<E on [r ,f). (A.4) 
Integrating (|A.4D we infer that there exist constants C+ > C_ > 0, such that 

C_(r - to) < t(r) - t < C+(r - r ) , \x l (r) - x | < SCI 1 ^ - l)' 1 ^ (A.5) 
on [t ,t). 

Making use of (|A.1() the geodesic equation (|A.3a|) yields 

d 2 t 

-4 = /,(r) where \l(t)\ < const p-*+*/ d , (A.6) 
dr z 

and by integration 

|^(r) - a-\r )\ < const gd+^O , hence |^ _ a -i| < const £<i+a/«*-fl ^ (AJ) 

where the constants depend on e, E, but are independent of to. The geodesic equation (|A.3b|) 
can be treated by noting that 

otrfc = (q - e dA ' 2 )tv~k + e dA ' 2 ^trfc + [3 + ^e - ^ 2 * -1 ) - [3 + ^] r 1 . (A.8) 

By (jA.lj) we obtain 

d 2 x? /_ 4\ ~_! d£ dx* 
dr 2 

integration leads to 



+ f 2 + -\ £-1^.^ = ? ( r ) W here | ? (r)| < const T 2 " 4 ^ ; (A. 
y dy dr dr 



& < const (t- 1 - 4 ^ - tot"- 2 " 4 ^ and f \ — \ < const ST 1 " 4 /** . (A.10) 
dr \ J dr 

If to is sufficiently large, then const t^ l+2 ^ d ® < e and const t 1 4 / d+< > < E, hence (|A.7j) 
and (lA~T0l) improve (|A.4|) on [to, t). Since r was chosen maximal, r must be infinite, and the 
above estimates hold globally. In particular, from (|A.10|) . 

1^1 < const t~ l - A l d and 1^ - cT 1 ! < const t~ 4 / d , (A.ll) 
dr dr 
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where the second inequality results from the fact that the geodesic is affinely parametrized. 
Global existence of the timelike geodesies orthogonal to the hypersurface t = to has thus been 
established; the asymptotic properties are given by (|A.11|) . 

To show that the family of geodesies gives rise to a global Gaussian coordinate system we 
investigate geodesic deviation. Consider the geodesic 7(7") and let dj 1 be the deviation vector 
between 7 and an infinitesimally neighboring geodesic in the direction dj, i.e., 

dd dd > 

d,°(r ) = 0, d/(r ) = 5 l j , -^-(to) = -(a~ 2 d ja )(r ) , -JL(t q ) = 0. (A.12) 

In analogy to l)A.4|) we can assume that {ddf/dr + a~ 2 dja\ < e and t^ddf/dr] < E on a 
maximal interval [tq,t). Along the lines of the above argument, by using derivatives of the 
geodesic equations (jA.3|) w.r.t. spatial variables we can improve these inequalities to obtain 
f = 00, and we get that 

dd i 

< const t~ x ~^ d so that \df - <5}| < const i~ A/d (A. 13) 

on [tq, 00). For large t the component d® will increase linearly in t. In line with the existing 
gauge freedom we may redefine d ■ 




for a suitable choice of \{t). We infer that the deviation vector behaves in a nice manner, at 
least for to sufficiently large, so that the family of geodesies originating from t = to generates 
a Gaussian coordinate system that is global in the future. ■ 

Let {t, x 1 } denote the Gaussian coordinate system constructed above. In these coordinates 
the field equations take the form d T (jij = —2kij, 

drtij = Kj + (trjfe)jb^ - 8tt^ + 47r5 i j (trS - p) , (A. 15a) 

-d 2 tp + A<p + (trk)d T tp = V'(tp), (A.15b) 

where Sij and p stem from the scalar field energy- momentum tensor (J2J), i.e., p = Too, 
Sij = Tij. 

It can be shown by example that series of the type 

W = T 2+A/d E $ V ~ am/d ' V = l^~ 2 (d + 2)/d]^ logr+E <P(m) r- 2m ' d (A.16) 

m>0 m>0 

do not provide solutions of (|A.15a|) in general. The asymptotic expansions of g, if, etc. 
necessarily include logarithmic terms, i.e., 

Lm 

m>0 /=0 

where L m G N Vm. It turns out that L m = for all m < d/2 in the case d = 4k, k G N, and 
L m = for all m < 3 + d in the case <i = 2(2A; + 1), fc € N. 
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It is interesting to contrast (|A.17|) and (JJJ) : in Gaussian coordinates the asymptotic expansions 
contain logarithmic terms in general, however, by the use of a suitable time coordinate these 
logarithmic terms can be removed. In Appendix[B]we investigate whether this statement can 
be generalized, on the level of formal power series, for arbitrary exponents d. 



B Formal asymptotic expansions with general exponents 

For a spacetime (M,g,{p) consider the Einstein scalar field equations with potential 



V(ip) = V exp ^-Kv^y *A with deR. (B.l) 

In [S] it was shown that the equations admit power series as formal solutions. The analysis 
was performed in Gaussian coordinates {r, x" 1 }, so that g = —dr 2 + gijdx l dx^, and it was 
found that 

W = S +A,d £ 9^~ 2m/d > V = l^~ 2 (d + 2)/^ 1/2 logr + £ <p (m) r" 2 "^ , (B.2) 

mGM mGAf 

where M is the set {0} U {{d/2)n\ + 2n 2 + (3 + d)n^ \ ri{ 6 N}. However, the series contain 
logarithmic terms, cf. (|A.17|) . if there exists n\,n% S N such that (d/2)m + 2n 2 = (3 + d), 
which is the case when d = n or d = 6/n, n £ N. 

The problem simplifies considerably when we make the ansatz 



~ g = e <$> (-dt 2 + g t3 dx l dxi) , if = k- 1 ^^ </> , (B.3) 

which is inspired by Let k % ■ denote the second fundamental form of the hypersurfaces 

t = const in the spacetime (M, —dt 2 + g), and o - * ■ its trace- free part. Then the Einstein scalar 
field (evolution) equations become dtg%j = — 2gnk and 



4a) 



= % + (trfc)a^- - i (v^Vrf - ivVV*^ - (Vv,> - ^V fc V fc <^ (B. 

<9 4 trfc = i? + (trA;) 2 - k; 2 V&-^— ^ - d^tik - ^V k (j)V k (j) - V k V k <t> (B.4b) 

□0 + V fe 0V fc - (d t ct>) 2 = -k 2 -^ V , (B.4c) 

where R and V refer to g. It can be proved that this system of equations, complemented by 
the constraints, admits power series of the following type as formal solutions: 

9ij = e 2Ht Uf + £ gff e~ mm ) , <P = dHt + (o) + £ (m) e""^ , (B.5) 

\ meM / m£M 

cf. (EU, where H 2 = 2V K 2 (d + 2)~ 1 (d + 3)" 1 and M = {2m x + (3 + d)m 2 \ rrn £ N}. The 
recursive algebraic system specifying the coefficients is essentially identical with (|35|) . the free 
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data is represented by the th and the (3 + d) th coefficients. In this context the exponent 
d € M is still not completely arbitrary, though: for d = 2k + 1, k 6 N, the expansions (|B.5|) 
must be supplemented by logarithmic terms. 

These results suggest that there exist two types of logarithmic terms in formal expansions: (i) 
artificial logarithms which are due to an unsuitable choice of coordinates and can be removed 
by using different coordinates, and (ii) genuine logarithms. In our particular case we have 
seen that, of the logarithmic terms in [S] which appear for d = n and d = 6/n, n € N, only 
those in the case d = 2fc + l, ft€N can be genuine. 
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